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ABSTRACT 

, We mapped 12 massive protostellar candidates in the CO 3=2—1 line, which in 

combination with Zhang et al. (2005) completes an unbiased survey of outflows for 
all 48 sources with />50° in a sample of 101 massive protostellar candidates. We 
detected outflows in 10 sources, implying 88% occurrence frequency of outflows for 
^ ■ the 48 sources. This supports the conclusion of previous studies that bipolar out- 

t^J- ■ flows are an integral component in the formation process of massive stars. The vast 

majority of the observed outflows are much more massive (>10 Mq) and energetic 
, (>100 Mq km s -1 ) than outflows from low- mass protostars. They also have large 

mass outflow rates (>2xl0~ 4 M Q yr _1 ), suggesting large (~lxl0~ 4 Mq yr _1 ) accre- 
^r) [ tion rates sufficient to overcome radiation pressure of the central massive protostars. 

We compared the frequency distribution of collimation factors of 40 massive outflows 
■ including those of this study with that of 36 low- mass outflows from the literature, and 

found no significant difference between the two. All these results are consistent with 
— . . the suggestion that massive stars form through accretion as do low-mass stars but with 

O . much higher accretion rates. 

o : 

Subject headings: ISM: clouds — ISM: jets and outflows - ISM: kinematics and dynamics 
— stars: formation 

s_i ■ 1. Introduction 

d; 

It has been well established that low-mass stars form by gravitational collapse of dense molecu- 
lar cores and that accretion disks and bipolar outflows are basic building blocks of the process (Shu 
et al. 1987). Compared to low-mass star formation, very little is known about the formation process 
of massive (>8 Mq) stars, such that it is still much debated whether massive stars form in a qual- 
itatively similar manner to low- mass stars. This debate is mainly because the radiation pressure 
of massive stars is strong enough to reverse materials infalling at the typical mass accretion rates 
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observed in low-mass young stellar objects (YSOs) (Wolfire & Cassinelli 1987). There are currently 
two major theoretical models in competition: accretion via disks (e.g., Yorke & Sonnhaltcr 2002; 
McKee &; Tan 2003) and coalescence of low-mass (proto)stars (e.g., Bonnell et al. 1998). The most 
convincing way to differentiate between the two models may be to investigate whether well-defined 
accretion disks exist around massive YSOs. But it is difficult to detect the accretion disks, if any, 
because they would be small (at most several hundred AU), short-lived, and easily confused by 
envelopes. In this study we attempt to distinguish between the two models by observing bipolar 
molecular outflows that are much easier to detect. 

Bipolar outflows from low-mass YSOs have been extensively studied since the early 1980's 
(see Bachiller 1996, and references therein), while systematic studies of outflows associated with 
massive star formation began much later. These studies are roughly divided into two groups. The 
first group are single-point CO line surveys towards samples of massive YSOs in search of high- 
velocity molecular gas (e.g, Shepherd & Churchwell 1996a; Sridharan et al. 2002). They showed that 
high-velocity gas is a common feature in massive YSOs. Shepherd & Churchwell (1996a) made a CO 
J=l— line survey towards 122 massive YSOs, the vast majority of which are ultracompact (UC) 
HII regions, and detected high-velocity gas in 90% of the 94 sources suitable for line wing analysis. 
Sridharan et al. (2002) obtained a similar detection rate (84%) from a CO J=2— 1 line survey of 69 
massive protostellar candidates. The second group are CO line mappings of subsets of the sources 
that exhibit high- velocity wings (e.g., Shepherd & Churchwell 1996b; Beuther et al. 2002a). These 
studies investigated whether the high-velocity wings are caused by bipolar outflows and derived the 
physical and dynamical parameters of the outflows. Shepherd & Churchwell (1996b) mapped 10 
sources with high- velocity gas and identified outflows in 5 sources. If we assume a similar detection 
rate for the remainder of their sample, the occurrence frequency of outflows would be about 45% 
for their sample. Beuther et al. (2002a) found outflows in 21 of 26 sources showing high-velocity 
wings in the sample of Sridharan et al. (2002), indicating a significantly higher occurrence rate 
(~70%). These studies also show that bipolar outflows from massive YSOs are much more massive 
and energetic than those observed in low-mass YSOs (see also Ridge & Moore 2001). There is some 
controversy as to whether massive outflows are much less collimated than low-mass outflows (see 
§3.1). 

The present paper presents CO J=2— 1 line maps of 12 massive star-forming regions at 27" 
angular resolution (Table 1). Combining the results of this study with those of Zhang et al. (2001, 
2005), we complete a census of molecular outflows for all sources with />50° in a flux- limited sample 
of massive protostellar candidates (see § 2.1), and estimate the occurrence frequency of molecular 
outflows for the sample. We also examine the morphological, physical, and dynamical properties 
of bipolar outflows observed, and compare them with those of outflows from low-mass YSOs. The 
selection criteria of the sources and the details of the observations are described in § 2. The results 
are presented in § 3. The implications of our results for massive star formation are discussed in § 4. 
We summarize the main results in § 5. 
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2. 



Observations 



2.1. 



Source Selection 



Our sources are a subset of the catalog of 101 massive protostellar candidates (Molinari et 
al. 1996), which were selected from a complete flux-limited sample of bright IRAS point sources 
with far-infrared (FIR) colors similar to compact molecular cores, based on the detection of NH3 
line emission. The original sample consists of 260 IRAS point sources that satisfy flux, color, and 
practical criteria (Palla et al. 1991): 1) F 60 >100 Jy and no upper limits for F25, F 60 , and Fi 00 , (2) 
0.61<log(F 60 /-F 25 )<1.74, and 0.087<log(Fi o/i ? 6o)<0.52, (3) no positional coincidence with known 
HII regions, and (4) |6|<10° and <5>-30°. Molinari et al. (1996) surveyed NH 3 (1,1) and (2,2) 
lines towards 163 of the 260 sources, and detected 101 sources. Among the 101 sources, 63 ("high" 
group) have colors characteristic of UC HII regions, log(.F25/-Fi2)>0.57 and log(i ? 25/i ? i2)^l-30 
(Wood & Churchwell 1989), while the remaining 38 ("low" group) have colors of log(F2s/Fi2)<0.57. 
Most (~90%) of the 101 sources have bolometric luminosities >10 3 L Q , so they are very likely 
to be massive protostars that have not yet formed UC HII regions, probably because of active 
accretion (Molinari et al. 1996, 1998). Zhang et al. (2001, 2005) mapped 69 of the 101 sources 
in the CO 3=2—1 line: 67 sources observed by Molinari et al. (1998) in radio continuum and 2 
other sources (IRAS 00117+6412 and IRAS 06056+2131). They detected outflows in about 90% 
(35/39) of their sources in the Galactic longitude range />50°. For sources at 0°</<50°, they could 
not compile reliable statistics because outflow signatures are frequently confused with other cloud 
components along the same lines of sight. In this study we observed 12 sources (Table 1). They 
include all 11 sources with />50° among the 34 sources not observed by Molinari et al. (1998), and 
IRAS 05274+3345. Thus our sample has 3 common sources (IRAS 00117+6412, IRAS 05274+3345, 
and IRAS 06056+2131) with the sample of Zhang et al. (2001, 2005). 



We mapped the 12 sources in the CO J=2— 1 line using the 12 m telescope at Kitt Peak 
in 2003 January and February. The telescope has a FWHM of 27" at 230 GHz. A 256-channel 
filterbank with 128 MHz bandwidth was used in parallel mode to observe both circular polarizations 
simultaneously, and the velocity resolution is 0.68 km s -1 . The system temperature varied in the 
range 500—700 K during the observing sessions. All the sources were mapped with 8" spacing in 
position switching mode using the on-the-fly observing technique. The reference positions were 
checked to be free from appreciable (>0.1 K) CO emission. They were usually ~30' away from 
the sources. The pointing and focus were checked every four hours with Venus, Saturn, or DR21. 
The typical rms noise level is 0.2 K. The observed temperature (T^) was converted to main-beam 



lr The 12m telescope is a lacility ol the National Science Foundation currently operated by the University of Arizona 
Steward Observatory under a loan agreement with the National Radio Astronomy Observatory. 



2.2. 



CO J=2— 1 Line Observations 
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brightness temperature (Tb) using the corrected main-beam efficiency (r?j^ b =0.51) provided by the 
observatory 

3. Results 
3.1. CO Maps 

We detected high-velocity CO gas towards all sources in our sample with the possible exception 
of IRAS 06105+1756 (Table 1). Here we compared the observed full width at zero intensity (FWZI) 
of each CO line with that expected from the FWHM by Gaussian fit, and define a high-velocity ex- 
cess if the former is greater than the latter. With the exception of one source (IRAS 22272+6358A) 
we find, as do Shepherd & Churchwell (1996a) in the CO J=l— line data, that the FWZI of the 
CO J=2— 1 line is always >15 km s _1 for the sources showing high-velocity excess. The FWZI of 
IRAS 22272+6358A is about 12 km s -1 in our 0.2 K rms data; possibly it would be broader at the 
0.02 K rms level of Shepherd & Churchwell (1996a). 

Based on contour maps of blue- and red-shifted high-velocity gas and position-velocity dia- 
grams, we conclude that 10 of the 11 sources host outflows; the exception is IRAS 06061+2151 (Fig. 
1). In IRAS 22272+6358A only blue-shifted high- velocity gas is detected. This is consistent with 
the result of CO J=l— line observations made by Sugitani et al. (1989) at 17" resolution. The 
IRAS sources are located reasonably close to the outflow centers in most cases, but they are offset 
from the centers by greater than the 27" telescope beamwidth in three sources (IRAS 06103+1523, 
IRAS 06382+0939, IRAS 23545+6508). In these cases the IRAS sources are unlikely to be the driv- 
ing sources of the outflows. Some maser sources and/or (sub)millimeter continuum sources, which 
are good indicators of massive star-forming sites in early evolutionary stages, are found around the 
outflow centers in IRAS 06382+0939 (Wolf-Chase et al. 2003) and IRAS 23545+6508 (Beuther et 
al. 2002b; see also Anglada & Rodriguez 2002). 

The morphologies of blue- and red-lobes are roughly symmetric in five sources (IRAS 00117+6412, 
IRAS 05274+3345, IRAS 06056+2131, IRAS 20227+4154, IRAS 22267+6244), whereas they are 
asymmetric in the remaining 4 sources (IRAS 05553+1631, IRAS 06103+1523, IRAS 06382+0939, 
IRAS 23545+6508). Such an asymmetry is often seen in the CO maps of outflows both from low- 
mass and high-mass YSOs (e.g., Bontemps et al. 1996; Beuther et al. 2002a; Zhang et al. 2005). 
Possible explanations of this feature are the asymmetric ejection of outflow-driving agents from 
the central stars, the nonuniform distribution of the ambient molecular gas, and the presence of 
multiple outflows. There seem to be two (or more) outflows in IRAS 06056+2131. This may be 
true for IRAS 05274+3345 as well (see also Hunter et al. 1995). This multiplicity complicates the 
high-velocity gas distributions of the two regions. We estimate the physical and dynamical param- 
eters of the two (A and B) outflows in IRAS 06056+2131 separately (see § 3.2), while we do not 
for the two of IRAS 05274+3345, because they are entangled in position as well as in velocity. 
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As noted earlier, there is some debate about the collimation of outflows from massive YSOs. 
Richer et al. (2000) and Ridge & Moore (2001) argued that outflows from massive YSOs are 
significantly less collimated than those observed in low-mass YSOs (see also Wu et al. 2004 and 
§ 4), while Beuther et al. (2002a) suggested, based on higher-resolution (11") CO maps, that they 
are better collimated than previously thought. We determine the collimation factors (/ c ), defined 
as the ratio of outflow length to width, of the 11 observed outflows (Table 2). The values range 
between 1.0 and 3.8. They are greater than 2.0 for 6 sources; for comparison, Ridge & Moore 
(2001) reported collimation factors of ~1— 2. Our average collimation factor is 2.3, similar to the 
2.1 value determined by Beuther et al. (2002a) for 15 outflows from massive protostellar candidates 
despite our factor of ~2.5 lower angular resolution. This may be because all outflows in our sample 
are located closer to the Sun and in less crowded (75°<Z<205°) regions of the Galactic plane than 
the majority of their sources. The average value is also very similar to that (2.2) measured for 36 
low-mass outflows as discussed in § 4. From our sample, IRAS 05274+3345 A, IRAS 05553+1631, 
and IRAS 06382+0939 are particularly well collimated with factors >3.0. In comparison, none of 
outflows observed by Beuther et al. (2002a) have collimation factors greater than 3.0. 



3.2. Outflow Parameters 



We derive the physical and dynamical parameters of the outflows in a similar way to Cabrit 
& Bertout (1992) and Beuther et al. (2002a). Here we assume (1) the excitation temperature of 
CO gas is 30 K (Molinari et al. 1996), (2) the integrated intensity ratio of CO to 13 CO is 10 in 
the line-wing velocity range (Choi et al. 1993), (3) [CO]/[H 2 ]=lxl0- 4 and [CO]/[ 13 CO]=89, and 
(4) the mean molecular weight is 2.3. We do not correct for the inclination of the outflow. Table 2 
displays the results: masses Mb (blue-shifted), M r (red-shifted), M out (total), momentum P, kinetic 
energy E, dynamical age i, mass outflow rate M ou t, mechanical force F m , mechanical luminosity 
L m , and collimation factor f c . The vast majority of outflows have sizes of 0.3—1.4 pc, dynamical 
ages of (1 — 10)xl0 4 yr, masses >10 M , momenta >100 M & km s -1 , and mass outflow rates of 
(1-10) xl0~ 4 M Q yr" 1 . Thus these outflows are much more massive and energetic than outflows 
from low-mass protostars (Class and Class I objects) with similar dynamical ages (e.g., Cabrit & 
Bertout 1992; Bontemps et al. 1996; Wu et al. 2004). This is consistent with the results of previous 
single-dish studies of massive YSO outflows (Shepherd & Churchwell 1996b; Beuther et al. 2002a; 
Zhang et al. 2005). The large estimated outflow rates suggest mass accretion rates large enough to 
overcome the strong radiation of the central massive stars (see § 4 for details). 



3.3. Occurrence Frequency of Outflows 

Zhang et al. (2001, 2005) mapped 69 of the 101 sources of Molinari et al. (1996) in the CO 
J=2— 1 line at about 30" grid spacing using the Kitt Peak 12 m and the Caltech Submillimeter 
Observatory 10 m. They detected outflows in 39 of the 65 sources that have data suitable for 
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outflow identification, for a detection rate of 60%. But it was practically impossible to compile 
reliable statistics for 26 sources in the Galactic longitude range 0°<Z<50° because of multiple cloud 
components in the same lines of sight. They found a much higher 90% (35/39) detection rate for 39 
sources with Z>50°. In this study we observed 12 sources selected from the same catalog. They are 
all located at />50°. We detected in 10 sources outflows that have similar physical and dynamical 
properties to other outflows from massive YSOs. This indicates the outflow detection rate is 88% 
for all 48 sources with />50° in the Molinari catalog, taking into account that the two samples have 
3 sources in common. It may be reasonable to expect a similar occurrence frequency of outflows for 
the sources at 0°<Z<50°, because both groups were selected by the same criteria. The FWZI of the 
CO line is ~10 km s _1 for the sources that show no distinct high-velocity wings. Assuming that 
the typical outflow velocity is ~25 km s _1 , outflows at inclinations >80° would not be detected. 
Statistically, the fraction of those outflows is <~10%. Thus this result strongly suggests that nearly 
all sources in the Molinari catalog have molecular outflows. 

There have been several attempts to measure the occurrence frequency of molecular outflows 
in low-mass protostellar objects (e.g., Terebey et al. 1989; Parker et al. 1991; Bontemps et al. 1996). 
The detection rates range between ~70% and 90%. For example, Bontemps et al. (1996) mapped 
45 low-mass protostars (36 Class I and 9 Class objects) of nearby (<450 pc) molecular clouds in 
the CO J=2— 1 line, and found outflows in 75%— 80% of the sources. The detection rate was higher 
(80%— 90%) for 34 YSOs with dust continuum emission. These values are similar to the detection 
rate determined by this study for the Molinari sources. Therefore, bipolar outflows appear to be 
an integral stage in the formation process of massive stars as well, as suggested by some previous 
studies (e.g., Zhang et al. 2001, 2005; Beuther et al. 2002a; Wu et al. 2004). 

4. Implications for Massive Star Formation 

As mentioned earlier, the formation process of massive stars is still poorly understood. Two 
major theoretical models are currently competing: accretion and coalescence. If massive stars form 
by accretion as do low-mass stars, all massive protostars would generate massive and powerful out- 
flows. Probably this is not true for the coalescence model where massive stars form from the collision 
of low-mass (proto)stars. There are three major collision units, such as stars, disks (/outflows), and 
cores, and six types of collisions are possible among the three (Bally &: Zinnecker 2005). Only if 
disk-disk collisions dominate the other types of collisions, and if the collisions do not profoundly 
disrupt the outflows of low-mass (proto)stars, might one expect that most massive protostars have 
massive and energetic outflows. Thus the observed high occurrence frequency of outflows can not 
be easily explained by the coalescence model. On the other hand, the collimation degree of outflows 
could be very different between the two models. The masses of massive outflows are usually one 
order of magnitude larger than those (<1 Mq) of low-mass outflows. If massive stars form by 
coalescence, this suggests that a massive outflow may arise after multiple collisions. In that case, 
although the outflows of low-mass protostars survive the collisions, the resulting massive outflow 



- 7- 



might be highly impulsive and poorly collimated like the OMC-1 outflow observed by Allen & 
Burton (1993). In contrast, one would expect relatively well-collimated outflows around massive 
protostars in the accretion model. 

Most outflows observed by this study do not seem to be very different from many low-mass 
outflows in collimation (e.g., Bontemps et al. 1996). In order to investigate this quantitatively, we 
compare the frequency distribution of collimation factors of 40 massive outflows with that of 36 
low-mass outflows. In this comparison we use the values derived by Lada (1985), Beuther et al. 
(2002a), and this study for 14, 15, and 11 massive outflows, respectively. The associated IRAS 
sources almost always have bolometric luminosities from 10 3 to 10 5 L , and the majority of them 
are massive protostellar candidates. For low-mass outflows we adopt the values determined by 
Lada (1985) for 12 sources, and estimate collimation factors for 24 sources using the CO maps 
provided by Myers et al. (1988), Schwartz et al. (1988), Parker et al. (1991), and Bontemps et 
al. (1996), which are systematic single-dish studies of multiple sources mainly at relatively low 
(30"— 45") resolutions. Here we exclude CO maps that do not cover the full outflows. Most of 
their driving sources are Class and Class I objects. Figure 2 shows the result. Though a higher 
fraction of low-mass outflows have collimation factors >3.5, the two distributions appear to be 
quite similar and the average values, 2.2, are the same. This does not agree with the finding of 
Wu et al. (2004) for a sample of 231 YSOs that the average collimation factors are 2.8±2.2 and 
2.1±1.0 for low-mass and massive outflows, respectively. This discrepancy may be mostly because 
of different resolutions of the data used for low-mass outflows in the analyses, even if we cannot 
exclude the possibility that some intrinsic difference exists between low-mass and massive outflow 
collimations. Our analysis includes only single-dish data for all low-mass outflows, while the Wu's 
analysis contains significantly higher-resolution single-dish and interferometer data for many low- 
mass outflows, which can reveal highly-collimated outflows. There are 10 outflows with collimation 
factors >5 in the Wu sample of low-mass outflows. Unless they are included, the average collimation 
factor is close to 2.2. On the other hand, it is not straightforward to identify highly-collimated 
outflows around massive YSOs even in interferometer data, because massive stars tend to form in 
groups and clusters where multiple outflows may overlap. For example, Beuther et al. (2002c) 
found at least 3 outflows, two of which are highly collimated, in the interferometer data of the 
massive star-forming region 05358+3543. The three could not be separated in single-dish data and 
the collimation factor of the combined outflow was measured to be 2.0 (Beuther et al. 2002a). 

It is widely accepted that low-mass outflows are momentum-driven by highly-collimated jets 
and/or wide-angle winds from the central (proto)stars (Konigl & Pudritz 2000; Shu et al. 2000). 
If massive outflows are produced by the same mechanisms, the mass accretion rates can be esti- 
mated from the observed outflow rates under the assumption of momentum conservation between 
jets/winds and outflows. The ratio of wind mass-loss rate to accretion rate is typically ~0.1 in the 
highly-collimated jet model (Konigl & Pudritz 2000) and ~0.3 in the wide-angle wind model (Shu 
et al. 2000). Assuming that the velocity ratio of winds to outflows is around 20 (see § 4.2 of Beuther 
et al. 2002a), the mass accretion rate would be 20%— 50% of the outflow rate. The average value of 
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the observed outflow rates is 5xlCT 4 M & yr~ 4 (Table 2), so the accretion rate is ~lxlCT 4 M & yr -1 . 
This value is much greater than the typical accretion rates of (10 -7 — 10 -5 ) M@ yr _1 expected and 
measured in low-mass YSOs (e.g., Shu 1977; Bontemps et al. 1996). It is also large enough to 
overcome the radiation pressure of the central stars, which have bolometric luminosities <2x 10 4 L & 
(Walmsley 1995). Therefore, the results of this study appear to support the accretion model rather 
than the coalescence model. 

5. Conclusions 

We mapped 12 massive protostellar candidates in the CO J=2— 1 line, which in combination 
with Zhang et al. (2005) completed a census of bipolar molecular outflows for all 48 sources with 
/>50° in the Molinari catalog of 101 massive protostellar candidates. We detected outflows in 
10 of the 12 sources, indicating 88% occurrence frequency of outflows in the 48 sources. The 
observed outflows are much more massive (>10 M@) and energetic (>100 M@ km s _1 ) than those 
observed in low- mass YSOs. These results imply that almost all sources of the Molinari catalog have 
very massive and powerful molecular outflows. Thus bipolar outflows appear to be a ubiquitous 
phenomenon in massive star formation as in low-mass star formation, as suggested by some previous 
studies. 

Nearly all outflows observed by this study have mass outflow rates >2xl0~ 4 M Q yr" 1 , sug- 
gesting large accretion rates of ~lxl0 -4 M@ yr -1 . These values are much greater than the typical 
accretion rates expected and measured in low-mass YSOs and are large enough to overcome the 
radiation pressure of the central massive protostars. We compared the collimation factors of 40 
massive outflows with those of 36 low-mass outflows, and found no significant difference between 
them. This can be naturally understood in the accretion model, but cannot be easily explained by 
the coalescence model. Therefore, it is likely that accretion plays a dominant role in the formation 
process of massive stars, although we cannot exclude the possibility that coalescence plays a role 
in some special circumstances. 

We acknowledge support from the Laboratory for Astronomical Imaging at the University of 
Illinois and NSF AST 0228953. 
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Table 1. Source Summary 
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a Source number and High/Low (H/L) designation of Molinari et al. (1996). 
b Wouterloot et al. (2004) and Brand et al. (1994). 
c Szymczak et al. (2000). 

d Assumed to be located in the Gemini OBI cloud complex (Carpenter et al. 1995), in which case the Lboi's were re-calculated with new distances 
from the values given by Molinari et al. (1996). 

° Assumed to be associated with the Cygnus X complex (Dame & Thaddeus 1985). 



Table 2. Outflow Properties Derived from CO J=2— 1 Line Data 
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a The properties of 05274+3345A and 05274+3345B are not separately estimated except for the collimation factor, because the two are entangled both 
in position and velocity. 

b Only blue-shifted high- velocity gas is detected (see the text). 



-13- 



00117+6412(2) 



i r 

[-50.0,-38.9] [-33.1,-20.1] 




0014 45 30 15 

RIGHT ASCENSION (J2000) 

Blue Contours: 3, 4, 6, 8, 10, 14, 18, 22 K km/s 
Red Contours: 3, 4, 6, 10, 14, 18, 22 K km/s 



05274+3345 (10) [-17.3, -6.2] [-0.4, 14.5] 




05 30 55 50 45 40 

RIGHT ASCENSION (J2000) 
Blue Contours: 6, 8, 10, 12, 15, 18, 21, 24, 27 K km/s 
Red Contours: 4, 6, 8, 10, 12, 14, 16, 20, 24, 28, 32, 36 K km/s 



05553+1631 (14) [-15.1,2.4] [8.9,25.2] 




05 58 25 20 15 

RIGHT ASCENSION (J2000) 
Blue Contours: 4, 6, 9,12,15,1 8, 21 , 24 K km/s 
Red Contours: 3, 4, 6, 8, 10, 12, 14 K km/s 



_ 06056+2131 (15) 


I I I 
[-14.9,-1.2] [6.6,20.2] 








A 






i 




~N i 


I I 



06 08 50 45 40 

RIGHT ASCENSION (J2000) 
Blue Contours: 4, 6, 8, 12, 16, 20, 24, 28, 32 K km/s 
Red Contours: 6. 8, 10, 12, 14, 16. 18 K km/s 



1 1 

06103+1523 (19) 



1 1 

[0.1,11.81 [19.6,35.2] 




06 13 25 20 15 

RIGHT ASCENSION (J2000) 
Blue Contours: 3, 4, 5, 6, 7, 8, 9, 10, 11 K km/s 
Red Contours: 4, 6, 9, 12, 15, 18, 21 K km/s 



"i i i i i r 
06382+0939(27) [-10.3,0.1] [10.5,20.3] 




06 41 15 10 05 00 40 55 

RIGHT ASCENSION (J2000) 
Blue Contours: 6, 8, 10, 12, 14, 16, 18, 22, 26 K km/s 
Red Contours: 6. 8, 10, 12, 14, 16, 18, 22, 26 K km/s 



Fig. 1. — CO J=2— 1 line integrated intensity maps. The Molinari source designation is given in parentheses 
following the IRAS name. Blue contours represent blue-shifted high-velocity gas, while red contours indicate red- 
shifted gas. The integrated velocity ranges are presented at the top right corner of each panel. Contour levels are 
listed at the bottom in each panel. The cross and its associated ellipse show the position of the IRAS point source 
and its uncertainty, respectively. The first and last panels show features marked by "N" which were not included 
when deriving the outflow parameters. In the first panel (IRAS 00117+6412) the "N" contours are a separate 
velocity component, whereas in the last panel (IRAS 23545+6508) the contribution of the "N" component is not 
well-determined because it is not fully imaged. 
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1. — Continued 
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Fig. 2. — Frequency distributions of collimation factors for (left) 36 low-mass outflows and (right) 
40 massive outflows. The values are listed in each bin. The two groups have the same average 
value, 2.2. 



